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THE  CURRENT  CONVECTIVE  INSTABILITY  AS  APPLIED 
TO  THE  AURORAL  IONOSPHERE 

I .  INTRODUCTION 

Recently  the  current  convective  instability  has  been  suggested  as  a 
possible  mechanism  responsible  for  irregularities  causing  scintillation 
phenomena  observed  by  the  DNA  Wideband  Satellite  during  diffuse  auroral  situa¬ 
tions  [Ossakow  and  Chaturvedi,  1979;  Chaturvedi  and  Ossakow,  1979;  Vickrey  et 
al . ,  1980].  The  observations  seem  to  suggest  the  scintillations  as  being 
caused  by  L-shell  aligned,  field  aligned,  large  scale  irregularities  co¬ 
located  in  the  regions  of  soft  particle  precipitation  and  horizontal  plasma 
density  gradients  with  the  source  region  confined  in  latitude  [Fremouw  et 
al . ,  1977;  Rino  et  al . ,  1978;  Rino  and  Owen,  1980].  The  current  convective 
instability  occurs  in  regions  where  the  field-aligned  currents  exceed  a  cer¬ 
tain  threshold,  determined  by  the  ambient  parameters,  etc.,  in  the  presence 
of  a  plasma  density  gradient  transverse  to  the  magnetic  field  [Kadomtsev  and 
Nedospasov,  1962]. 

The  linear  theory  of  the  instability  as  applied  to  the  auroral  situation 
appears  to  give  plausible  growth  times,  regions  of  occurrence,  etc.  [Ossakow 
and  Chaturvedi,  1979;  Vickrey  et  al ■ ,  1980],  Some  of  the  observed  features, 
like  L-shell  alignment  and  percentage  fluctuation  of  irregularity  amplitudes, 
may  be  explained  by  a  nonlinear  saturation  mechanism  of  the  instability  that 
invokes  mode  coupling  to  damped  modes  in  a  two-dimensional  treatment  perpen¬ 
dicular  to  the  magnetic  field  [Chaturvedi  and  Ossakow,  1979].  The  field- 
aligned  currents  cause  weak  shear  in  the  magnetic  field  with  shear  lengths 
comparable  in  order  of  magnitude  to  the  parallel  scale  sizes  of  the 
instability-driven  unstable  modes.  The  effect  of  shear  is  to  reduce  the 
growth  rates  somewhat  but  to  spatially  localize  the  mode  in  the  density 
gradient  direction  [Huba  and  Ossakow,  1980].  The  first  global  numerical 

simulations  [Kesklnen  et  al . ,  1980]  performed  for  the  nonlinear  current 
Manuscript  submitted  October  24,  1980. 
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convective  instability  show  that  during  the  nonlinear  evolution  the  regions 

of  enhancements  move  equatorward  while  those  of  depletions  drift  polewards. 

In  addition,  the  power  spectrum  of  the  density  fluctuations  follows  a  power 

_2 

law  with  a  typical  value  in  the  vicinity  of  ~  k 

In  the  present  work  we  have  extended  the  previous  linear  analysis  of  the 
current  convective  instability  [Ossakow  and  Chaturvedi,  1979]  to  include  ion 
inertial  effects,  which  assume  importance  at  high  altitudes  where  the  ion- 
neutral  collision  frequency  gets  smaller.  We  find  the  growth  rates  to 
remain  of  the  same  order  of  magnitude  in  the  inertial  case  as  they  were  in 
the  collisional  case  [Ossakow  and  Chaturvedi,  1979].  At  the  same  time,  we 
present  the  growth  rate  expression  of  the  current  convective  Instability  for 
the  other  limit  when  ions  are  strongly  collisional,  i.e.,  v^n  »  as  is 
the  case  in  the  E-region  ionosphere.  Further,  as  is  well  known,  in  consider¬ 
ing  long  wavelengths  such  that  ck/u)^  ^  1  (where  c  is  the  velocity  of  light, 
k  the  wave  number  and  u)g  the  electron  plasma  frequency),  one  should  consider 
electromagnetic  effects  also  [Kaw  et  al . ,  1974].  We  find  these  effects  to  be 
very  small  in  the  ionospheric  situation. 

In  the  next  sections,  we  first  present  the  growth  rate  calculations  for 
the  ion-inertia  dominated  regime  and  for  the  case  of  strongly  collisional 
ions,  then  the  effect  of  electromagnetic  corrections  on  the  growth  rate,  and 
in  the  final  section  we  discuss  them  for  the  ionospheric  application. 

II .  THEORY 

A.  Ion  Inertia  Effects 

We  shall  work  in  a  coordinate  system  which  has  its  z-axis  aligned  with 
the  magnetic  field,  1^,  with  its  y-axis  pointing  northwards,  and  the  x-axis 
pointing  westward.  A  field-aligned  current,  along  the  z-axis  and  a  plasma 
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density  gradient  along  ttie  y-axis  are  assumed  as  a  part  of  the  equilibrium. 
This  neglects  a  weak  altitude  dependent  gradient,  any  east-west  gradient  and 
any  equilibrium  electric  fields.  The  temperature  effects  responsible  for 
diffusive  damping  are  ignored  for  simplicity  as  is  magnetic  field  shear. 

The  electron  inertia  is  neglected  for  the  low  frequencies  concerned  in  the 
paper  and  ven»  the  electron-neutral  collision  frequency,  is  neglected  in 
comparison  to  v  the  electron-ion  Coulomb  collision  frequency,  which  is  a 
valid  assumption  at  F-region  altitudes.  The  zero-order  current  is  assumed  to 
be  caused  by  an  equilibrium  drift  of  electrons  over  ions  along  the  z-axis. 
The  perturbations  are  assumed  to  vary  as  “  exp  (ik^x  +  ik^y  +  ik zz  -  iivt)  and 
the  local  approximation  is  used  in  the  analysis  which  is  justified  as  long  as 
we  consider  wavelengths  to  be  small  compared  to  the  density  gradient  scale- 
length  (ky  »  d  In  nQ/dy)  . 

With  these  assumptions  the  electron  velocities  are  given  as 


eE 

v  -  -  —5- 

ze  mv 
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cE^xz 
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— Xe  B 
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(2) 


where  the  electrostatic  assumptions  for  perturbed  electric  fields, 
E  -  V <p  has  been  used. 

The  ion  momentum  transfer  equation  is 
d  v  , 

dt  -x  c  -xi  -ox  in  — xi 


(3) 


The  perturbed  ion-velocities  in  the  perturbed  fields  can  be  written  as 


^xi 


(A) 
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The  difference  of  eqs.  (8)  and  (9)  leads  to 
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where  £1^  »  v  has  been  assumed  in  addition  to  the  quasi-neutrality 
assumption  ng  <=»  n^  =  n. 

The  dispersion  relation  is  obtained  from  (8)  and  (10)  and  is 
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where  some  small  terms  have  been  neglected.  Splitting  uu  »  +  iy,  the 

growth  rate  is  given  as 


where 


(12) 


Equation  (12)  can  be  maximized  with  respect  to  k^/k^  and  this  maximum  growth 
rate  is  obtained  from  the  cubic  equation  for  k^/k^ 


AzY  /kzVeiVin  1  f^ei! 
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0. 


(12a) 


It  is  simple  to  recover  the  collisional  growth  rate  from  (12),  for 
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(13) 
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which  is  the  same  as  before  [Ossakow  and  Chaturvedi,  1979;  Chaturvedi  and 
Ossakow,  1979].  The  maximum  growth  rate  in  this  case  occurs  for  (which  also 
results  from  (12a)  by  neglecting  the  last  term  on  the  lefthand  side  of  (12a)) 
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In  the  limit  ="  0,  we  recover  the  purely  inertial  case  from  the 
growth  rate  expression  given  by  (12), 


.  k  2  0  n.  .  fk  4  Q  *0  2  4v  eko  ' 
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A  simplified  expression  of  growth  rate  for  this  case  can  be  obtained  when 
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then  one  gets 
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Maximum  growth  rate  in  this  case  occurs  for  (from  (15)  and  also  results  from 
(12a)  by  neglecting  the  second  term  on  the  lefthand  side  of  (12a)) 
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(16) 


It  should  be  noted  that  in  the  inertial  case  the  value  of  k  /k  which  maxi- 

z  x 

mizes  the  growth  rate  depends  on  instability  driving  parameters,  vq  and  e, 
and  not  just  background  Ionospheric  parameters,  Vg^,  v^,  ftg,  and  0^,  as  well 
the  case  in  the  collision  dominated  regime. 


B.  Coliisional  Ions 

We  now  consider  a  situation  wherein  ions  are  highly  coliisional  such 
that  »  ft  .  This  situation  is  frequently  realized  at  E-region  alti¬ 
tudes  in  ionosphere.  The  ion  and  electron  velocities  may  be  written  as 
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Upon  using  the  electrostatic  assumption,  one  can  readily  write  the  ion  and 
electron  continuity  equations  as 
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where  the  quasi -neutral ity  assumption  has  been  used.  The  difference  of  the 
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above  two  equations  yields 
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The  dispersion  relation  can  be  obtained  from  eqs.  (19)  and  (21), 
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One  puts  ti)  =  +  iy,  then 
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(22a) 


The  growth  rate  expression  from  eq.  (22)  is 
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(22b) 


For  typical  E-region  situations,  m/M  ~  2  x  10-5,  v  ~  2.5  x  103s  *, 

’  in 

Vgn  ~  4  x  lO^s-1,  and  assuming  L~10kms,vo  ~  2xl0Vm/s,  etc.,  we  find 

V ,  2  /V 

—  --^-n  r  Z4  »( rV  I  »  so  that  uj_  ~  3.5  x  10-2s-1  and  y  ~  10-3sec-1  ,  where 
mvk2\ll.  kL/  R 

en  x  \  i  x  / 

k  /k  ~  (M/m  — —  |  ^  ~  10-2  was  used  with  A  ~  300  m.  We  note  here  that  at 

Z  \  V  /  X 

\  en  / 

E-region  altitudes,  an  ambient  electric  field  transverse  to  the  magnetic 
field  causes  growth  of  perturbations  due  to  the  so-called  E  x  B  instability 
The  growth  rate  expression  for  this  cross-f ield  instability  may  be  written  as 
(Sudan  et  al . ,  1973) 

/k  e  n 
j! _  /  _? 

(1 
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+  Iv-en'V 
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VenVin  1  ^n0  cEq 

where  u>R  -  k-y^/d+4.)  ,  *  ~  ~  iT  df  ’  V"  =  B~ 

e  1  o  o 

(with  bein'  he  primary  transverse  electric  field).  Here  temperature 
effects  have  been  ignored  for  simplicity.  The  growth  rate  due  to  the  It  x  b 

instability  can  be  much  higher  depending  upon  the  parameters.  Thus,  for 

-1  -i  -3  -1 

v  ~  100  m/s,  L  ~  10  km,  y„  „  ~  10  s  ,  in  contrast  to  y  ~  10  s 

O  CjXd  CC 

obtained  above  for  the  current  convective  instability. 


C.  Electromagnetic  Effects 

In  this  section  we  present  results  on  the  electromagnetic  effects  on  the 
growth  rate  of  the  current  convective  instability.  We  consider  the  case  in 
which  temperature  effects  are  ignored  in  addition  to  the  ion  and  electron 
inertia.  Instead  of  the  electrostatic  approximation,  we  use  the  full  set  of 
Maxwell's  equations. 


.SB  ,  .  SE 

VxE  =  -  i  ;  VxB  =  —  J  +  -  ^r 
c  at  —  c  —  c  at 


The  two  equations  above  can  be  combined  into  one  as, 
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where  the  displacement  current  has  been  ignored  for  the  low  frequencies 
concerned.  The  perturbed  current  is 
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where  the  perturbed  ion  and  electron  fluid  velocities  are  given  by 
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The  rest  of  the  perturbation  analysis  follows  the  procedure  outlined  in 
the  preceding  section.  Use  of  (25)  and  (26)  in  (24)  gives  .J  which  then  is 
substituted  back  in  eq.  (23)  to  yield  three  component  equations.  The  dis¬ 
persion  relation  is  obtained  from  the  determinant  of  these  three  component 
equations  and  is,  after  some  algebra, 
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The  electrostatic  limit  can  readily  be  obtained  from  eq.  (27)  by  using  the 

U)  U) 

e  6 

approximation  — jj—  «  1.  In  the  opposite  limit,  1,  electromagnetic 

effects  can  be  stabilizing  or  destabilizing  depending  on  the  parameters. 
However,  for  the  auroral  F-region  situation  they  turn  out  to  be  negligibly 
small.  The  growth  rate  in  this  regime  is,  from  eq.  (27), 


(28) 


For  typical  parameters,  u)g  ~  1.7xl07  rad  s  1 ,  X^  ~  1  km,  u^/ck^  ~  10, 

v,  /Q.  ~  2.5x10  v  /Q  ~  10  4 ,  Q  ~  5xl06,  v  ~  SxlO^cm/s,  X  ~  lO^km, 
in  i  e  e  e  o  z 

L  ~  50  km,  the  second  term  in  eq.  (28),  though  stabilizing,  is  ~  10  6s  1 

compared  with  the  typical  growth  rates  of  the  current  convective  instability 

”3  ""l 

in  the  auroral  ionosphere  which  go  as  ~  10  s  .  Thus  we  find  that  in  the 
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ionosphere  the  growth  rate  is  unaffected  by  electromagnetic  effects  for  the 


current  convective 
the  electromagnetic 
destabilizing,  with 

Y  =“ 


instability.  In  a  parameter  regime  where  v  /Q  > 
effects  on  the  current  convective  instability  can 
the  growth  rate  becoming. 


be 


(28a) 


This  parameter  regime  may  be  satisfied  at  lower  altitudes  in  the  ionosphere. 


III.  DISCUSSION 

We  have  presented  in  the  previous  sections  the  growtli  rate  expressions 
in  the  case  when  ion-inertia  is  important  (eqs.  (12),  (15),  (16));  ions  are 
highly  collisional  (eq .  22)  and  in  the  case  when  the  electromagnetic  effects 
are  included  (eq.  (28)).  This  extends  our  previous  analysis  wherein  we 
applied  the  current  convective  instability  to  the  diffuse  auroral  situation 
ignoring  ion-inertia, and  using  the  electrostatic  approximation  in  the 
analysis  (eq.  (13),  ( 1 A ) )  and  [Ossakow  and  Chaturvedi,  1979]. 

For  typical  ionospheric  K- region  parameters  at  auroral  latitudes  at 

350-400  km  altitudes,  we  have  vq  ~  500  m/s,  m/M  ~  3x10  ',  v  ^  ~  5x1 0‘ , 

v.  ~  5x10  ,  c  1  ~  I.  ~  50  km,  we  see  that  from  eq .  (14),  v  ~  3x10  'sec  1 

in 

and  similarly  from  eq .  (16),  ~  3x10  'sec  *.  We  find  that  the  growth  rate 

remains  the  same  order  of  magnitude  in  the  collisional  and  inertial 
domains.  At  higher  altitudes,  where  Vj  10  \s  v  ^  ~  30,  one  finds  that  the 
growth  rate  for  the  same  values  of  v^  and  r  ',  is  ~  3x10  ^s  Thus  one  can 
conclude  that  the  current  convective  instability  is  applicable  equally  well 
to  higher  altitudes  in  the  ionosphere  where  becomes  smaller  in  contrast  to 


•<o< 


the  collisional  case  presented  in  our  previous  work  [Ossakow  and  Chaturvedi, 

1979].  Thermal  diffusion  will  damp  out  the  modes  at  a  rate  proportional  to 
Ve  kz2 

~  cs2^j.2  +  ^ —  )  [Chaturvedi  and  Ossakow,  19793  which  determines  the 

e  i  in 

threshold  for  the  instability,  and  for  a  given  set  of  parameters  correspond¬ 
ing  to  instability,  determines  the  cut-off  scale  lengths,  below  which  the 
system  is  stable.  For  the  present  set  of  parameters,  the  stable  (damped 
modes)  wavelengths  correspond  to  X^D  ^  70  m  and  XjjD  ^  400  km.  In  the  limit 
of  highly  collisional  ions,  corresponding  to  E-region  altitudes,  the  growth 
rate  can  be  on  the  order  of  ~  10  sec  .  However,  the  presence  of  transverse 
electric  fields  of  moderate  strengths  (a  few  mv/m)  would  cause  growth  at  a 
much  faster  rate  due  to  the  well-known  ExB  instability.  The  electromagnetic 
effects  on  the  growth  rate  can  be  similarly  estimated,  and  for  the  above  set 

-3 

of  parameters,  for  wavelengths  ~  1  km,  they  turn  out  to  be  a  factor  ~  10 
smaller  compared  to  the  electrostatic  growth  rate,  and  are  rather  small  to 
have  any  effect. 


This  work  was  supported  by  DNA  and  ONR. 
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